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A Method for the Reduction of Comet Photographs» 

By H. C. Plummer, M.A. 

1. The magnificent results which have lately been obtained in 
comet photography suggest some reflections on the method by 
which the plates may be reduced. The method of Bessel,* 
which was designed for the treatment of Halley’s comet in 1835, 
and has been followed more or less faithfully by all subsequent 
writers, deals with observed position-angles and distances with 
respect to the nucleus. If the observations are not made in this 
form, they must first be reduced to it. Of course position-angles 
and distances can be very easily deduced from the photographic 
plate, f But it has frequently been found | that there is an 
advantage in the use of rectangular coordinates over spherical 
geometry in dealing with photographs. There is a corresponding 
change in the method of calculation from the general use of 
logarithms to the calculating machine which appeals strongly to 
some computers. It may, therefore, be worth while to investigate 
the formulae which seem to be the natural expression of this point 
of view. Whether the emanations from a comet remain strictly in 
the plane of the orbit is open to doubt, but provisionally it is 
necessary to make this assumption, and it forms the basis of what 
follows. 

2. It will be convenient to give precepts for certain preliminary 
operations. Let us form the arrays 

— sin A cos A o 

- sin D cos A — sin D sin A cos D . (1) 

cos D cos A cos D sin A sin D 


where A, D are the R.A. and declination of the plate-centre, and 


X 

a sin(a + u ) 
acos(a + u) 


Y 

b sin Q 3 + u) 
b cos ( f 3 4- u) 


Z 

c sin (y + u) 
c cos (y + u ) 


<*> 


where X, Y, Z are the rectangular coordinates of the Sun referred 
to the equatorial system of coordinates, a, b , c, a, /3, y the Gaussian 
constants for the orbit of the comet, and u an angle measured from 
the ascending node in the direction of the comet’s motion: the 
actual value to be assigned to u can be left for definition later on. 
Next we form the array 


h 

h 

h 


m l 

n i 





m 2 

n 2 

. 

. 

• 

. 

m B 

n 3 






(3) 


* A.K, No. 300. t Of. M'N., lxix. p. 100. 

t Cf. M.N., lx. pp. 176 and 201 ; lxii. p. 29. 
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which is formed by the products by rows from the arrays ( 1 ) and 
( 2 ), so that 

l Y = - sin A . X + cos A . Y 
m 1 = - sin A. a sin (a 4- u) 4 - cos A . b sin (fi + u) 
n l = - sin A. a cos (a 4- u) + cos A. b cos (ft 4- u) 
l 2 — - sin D cos A. X - sin D sin A. Y 4- cos D . Z 
m 2 = - sin D cos A . a sin (a + u) — sin D sin A. b sin (/3 4 - u) + cos D . c sin (y 4 - u) 
n 2 = - sin D cos A . a cos (a 4 - u) - sin D sin A. b cos (/3 4 - u) 4- cos D . c cos (y 4 - u) 
cos D cos A. X + cos D sin A. Y + sin D . Z 

cos D cos A. a sin (a 4- u) 4- cos D sin A. b sin (fi 4 - u) 4 - sin, D . c sin (y 4- u) 
cos D cos A . a cos (a 4- u) 4- cos D sin A . b cos ([3 + u) 4- sin D . c cos (y 4 - u) 


n 3 = 


Lastly we form the conjugate array 


Li 




m 2 

n 2 

L s 

m 3 



(4) 


the constituents of which are the minors 


L 1 = m 2 rc 8 - m s M 2 , M 1 = l 2 n 3 -l 3 n 2 , N 1 = l i m s -l s m 2 
L- 2 = , M 2 -=l^n z - Z 3 w 1 , X 2 = Z 1 m 3 - 

L 3 — ~ > M 3 = l\U 2 — l 2 n x , X 3 — l 1 m 2 — 

The operations involved in computing the constituents of the 
arrays ( 3 ) and ( 4 ) will be very easily carried out with the help of 
an arithmometer. 

3 . Let p , 0 be the polar coordinates of a point in the plane 
of the orbit, the origin being the San and the initial line making 
an angle u with the direction of the node. Then the geocentric 
coordinates of the point are 

x j =-• ap sin (a 4 - u + 0) 4- X 
*= ax sin (a 4 - u) 4- ay cos (a 4- u) 4- X 
y l = bx sin (f3 4- u) 4- by cos (/3 4- u) 4- Y 
z x = cx sin (y 4- u) 4- cy cos (y + u) 4- Z 

But if x 2 , y 2 , z 2 are the coordinates of the same point referred to 
the axes 0 A, OB, 0 C, where C is the plate-centre, and the plane 
0 C B passes through the pole P, 

x 2 = — sin A . x 1 4 * cos A . y Y 

y 2 = — sin D cos A.. x x - sin D sin A . y x 4 - cos D. z 1 

z 2 = cos D cos A. x 1 4 - cos D sin A. y 1 4 - sin D . z l 

the coefficients being given by the array ( 1 ). Hence 

x 2 = l l + m x x 4 - n x y 
y 2 = h + m 2 x + n 2 V 
Z 2 ~ h + m 3 X d- n Zy 

But if i, rj are the standard coordinates of the image of the point 
on the plate, / is the focal length, and parallax is neglected, 

?/»2 = ’?/ y 2=//*2 


© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 


Downloaded from http://mnras.oxfordjournals.org/ at Dalhousie University on May 17, 2015 






1909MNRAS..69..191P 


Jan. 1909. Reduction of Comet Photographs. 193 

so that x , y , the rectangular coordinates in the plane of the orbit, 
are connected with £, y by the relations 

0 h + + w i2/)/f = (h + = {l 3 + m i x+ n s y)/f (5) 

which are equivalent to 

Ljf- L 2 , + L 3 /= (- M,£ + M a iy - M s f)/x = (Njf - N s , + N s f)/y (6) 

Thus the sets of coefficients in (3) or (4) enable us to pass easily 
from points in the plane of the orbit to the corresponding points 
on the photograph, or vice versa. 

4. The value to be assigned to u has purposely been left 
arbitrary. If, for example, we put u = o, the axis of x will be the 
line from the Sun to the ascending node of the comet. But if we 
wish to form a representation of a number of different points 
connected with the comet at the date of a particular photograph, 
it is natural to make the radius vector of the nucleus the axis of x , 
and then we must put 

u = vs — O + v 

where v is the true anomaly of the nucleus at the date of the 
photograph. On the other hand, if we wish to follow the actual 
motion of a particular point as shown by a series of photographs at 
different times, it will be advantageous to use a fixed axis of x by 
putting 

u = VS - O + v Q 

where v 0 is a constant chosen for convenience, and may be the true 
anomaly of the nucleus at the time of one of the photographs, or 
may be chosen for some particular reason; e.g. if we put v 0 = o, 
the axis of x will be drawn from the Sun in the direction of 
perihelion. In any case the Sun is at the origin of coordinates. 

5. The projective relation between corresponding points on the 
photograph and on the plane of the comet’s orbit make it easy to 
transform any curved aggregate of points on the photograph into 
the corresponding curve in space by means of (5). Thus if we 
observe a straight streamer whose equation on the plate is 

A£ + B?7 + C = o 

this corresponds in the plane of the orbit to 

+ B^ + CZg + + BW&2 + C77ig)fl; + "1“ = o 

And similarly any more complicated curve on the plate for which 
the equation can be found can be transferred to the plane of the 
orbit with far greater ease than by computing a series of positions 
point by point. 

6. When a particular condensation in the matter emitted by a 
comet can be observed on a series of plates, it may be possible to 
determine the motion, and hence to infer the law of force which is 
acting. The difficulties in the way of identifying such a conden- 


© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 


Downloaded from http://mnras.oxfordjournals.org/ at Dalhousie University on May 17, 2015 




1909MNRAS..69..191P 


194 Method for the Reduction of Comet Photographs, lxix . 3 , 

sation and measuring it with any approach to accuracy are very 
great. Hitherto attempts to do this have not met with very great 
success. But it is a problem of immense interest. If the equation 
of the orbit is 

£> = />{ 1 +e cos (0 -y )} 

and we have a series of positions (x r , y r ) at times t r , of which the 
number must be at least three, the constants will be given by the 
equations 

p = p r + ex r cos y + ey r sin y 

for these will determine p, e cos y and e sin y. When these have 
been found, three cases may arise according as (1) e<i, (2) e>i 
and p is positive, (3) e>i and p is negative, these corresponding 
to motion on an ellipse, on the concave and on the convex branch 
of an hyperbola respectively. 

7. It may be convenient to set down here the formulae which 
suffice to determine the constant of force in each of the three cases. 

(1) x = pcoaO, y — psin 0 , v — 0 -y 
e = sin cf> , a =p sec 2 </> 

tan -|-E = tan (45 0 - -|</>) tan \v 
p)a~§(t - T) = E - e sin E 
/x = k 2 (i-/x') 

which determines p, the ratio of the repulsive action to gravita¬ 
tion, 7 c being the Gaussian constant. 

(2) x = p cos 0, y =p sin 0, v = 0 -y 
e = sec <jf>, a = p cot 2 4 > 

tan JF = tan \cj> tan \v 

* p$a-%{t - T) = e tan F - Log tan (45° + ^^ 

p = k\ i-fi') 

which determines p! as before. Log denotes natural logarithm. 

(3) For clearness, when p as determined in §6 is negative, we 
simply change its sign and rewrite the equation of the convex 
branch 

p — - p + ep cos (0 - y — 7r) 

and then 

x = p cos 0 , y = p sin 6 , v = Q - y - ir 

e = sec <f ), a = p cot 2 <£ 

tan -|-G = cot tan \v 

p$a~l{t - T) = Log tan (45 0 + JG) + e tan G 

which again determines p! in terms of gravitation. Astronomical 
units are assumed throughout. 

8. If three observations of a condensation are available, it is 
evident that only a pair of times is required to determine p and T, 
and that the time of the third observation serves as a control. A 
more satisfactory control must be looked for in a greater number 
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of observations, and the success with which they can be brought 
into harmony. If discrepancies present themselves, there are only 
too many causes which, may contribute to produce them. These 
include not only false identifications and want of accuracy in a 
particularly difficult class of measurement, but also the possibly 
false assumption that the matter is moving in the plane of the 
orbit, and an error in the hypothesis that the repulsive force varies 
as the inverse square of the distance. It is to be hoped that 
those who have reaped such a brilliant success in the field of 
photographic observation will be undaunted by the patent 
difficulties, and will arrive at some more definite conclusions as 
to the laws of cometary forms than we can be said to have reached 
at the present time. 

University Observatory , Oxford: 

1909 January 7. 


The Clouds of Venus and their Significance. By Arthur W. 
Clayden, M.A., Principal of the .Royal Albert Memorial 
University College, Exeter. 

If it be assumed that the visible surface of Yenus is the outer 
face of a heavily cloud-laden atmosphere, comparable in mass and 
composition with our own, it seems that certain deductions may be 
based on its observed features which have not hitherto been pointed 
out, though they convey valuable hints as to the planet’s physical 
condition, the position of its axis, and the period of rotation. 

There is no need to recapitulate here the many reasons for be¬ 
lieving that the surface visible to most telescopes is actually cloud, 
and the sequel will show that such a view is not incompatible with 
the detection of sharp markings such as have been described by 
specially favoured observers. My own observations have been 
made with the 6*8 inch refractor by Tulley, described in vol. ii. of 
the Society’s Memoirs. Its long focus (144 inches) gives excellent 
definition of Jupiter and Saturn, when the seeing is good, with 
powers up to 450 ; but it has never shown me anything on Yenus 
which does not seem easily explicable on the assumption that the 
planet is shrouded in clouds, and that its atmosphere bears some¬ 
thing like the same relation to its mass as our own atmosphere 
does to the Earth’s. 

Taking this for granted, let us consider how those clouds should 
be distributed according to different views as to the period of rota¬ 
tion and the position of the polar axis. 

First let us assume that the axis makes a large angle with the 
plane of the orbit, and that the period of rotation is coincident with 
the revolution. That is to say, that the planet always turns the 
same face t© the Sun—one side always exposed to the solar glare, 

14 
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